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The quest for gravitational waves from coalescing binaries is customarily performed by the LIGO- 
Virgo collaboration via matched filtering, which requires a detailed knowledge of the signal. Com- 
plete analytical coalescence waveforms are currently available only for the non-precessing binary sys- 
tems. In this paper we introduce complete phenomenological waveforms for the dominant quadrupo- 
lar mode of generically spinning systems. These waveforms are constructed by bridging the gap 
between the analytically known inspiral phase, described by spin Taylor (T4) approximants in the 
restricted waveform approximation, and the ring-down phase through a phenomenological interme- 
diate phase, calibrated by comparison with specific, numerically generated waveforms, describing 
equal mass systems with dimension-less spin magnitudes equal to 0.6. The overlap integral between 
numerical and phenomenological waveforms ranges between 0.95 and 0.99. 



O I. INTRODUCTION 
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>>■ The Laser Interferometer Gravitational-wave Observatory (LIGO) and Virgo constitute a network of kilometer-scale 
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interferometers for the detection of gravitational waves. The initial detector configuration has successfully acquired 
data at design sensitivity; substantial upgrades over the next few years, to Advanced LIGO and Advanced Virgo, 
*K^ are expected to yield the detection of gravitational waves from the coalescence of compact binary systems of neutron 

stars or black holes. 

The coalescence of binary systems is usually described in terms of three distinct phases: the inspiral, the merger 

and the ring-down. The inspiral phase allows for an accurate analytical description via the so-called Post-Newtonian 

I (PN) expansion, see for instance [ij for a review. The ring-down also admits a perturbative analytical model, as it 

tp, describes the damped oscillations of the single object resulting from the binary coalescence, as a superposition of 

^_i black- hole quasi normal modes [2 . The merger phase is however fully non-perturbative and for generic systems it has 

not been described analytically but rather by numerical simulations. During the last five years numerical relativity 

^^ has made tremendous progress in describing the full coalescence of a binary system beginning with [3]-[5] and more 

recently [Blll2|. see |13H15) for reviews and it can now produce waveforms for generic spin orientations, with moderate 

spin magnitude (< 0.9) and mass ratios (< 10 : 1). 

\f\ Matched filtering is typically used in LIGO- Virgo data analysis in order to uncover weak signals buried into noise 

f^ [TB] . This technique requires a large number (tens of thousands) of templates to be checked against real data. Due to 

I _* the computational cost of numerical simulations, this is only possible with an analytical knowledge of the templates. 

(^— ^ Analytical methods are currently available to reproduce the complete waveform emitted by non-precessing coalescing 

f^ binaries. This has been achieved in the Effective One Body (EOB) construction pTHlQ] for non-spinning systems, 

T— I and in the EOB-spin waveforms ^^, an extension of the EOB method to the case of binaries with non- vanishing 

L[ spins aligned with the orbital angular momentum. Both spinning and non spinning EOB require a comparison 

• ^H with numerically generated waveforms in order to calibrate some free parameters of the model. Another method for 

^^ generating analytical waveforms for spinning non-precessing binaries is by joining PN-generated inspiral and numerical 

J-^ waveforms |10l I21j to construct phenomenological waveforms. 

In this paper we present a new family of analytical waveforms from the coalescence of generic spinning binaries. 
The waveforms are obtained by interpolating between the perturbative PN inspiral description and the ring-down, 
both admitting physically motivated analytic models, with a phenomenological phase over the merger portion of 
the signal. Since in the matched-filtering search an accurate determination of the gravitational phase of the signal is 
crucial, more important than an accurate amplitude determination, we focus on the phenomenological parametrization 
of the gravitational wave phase, and give lower priority to the amplitude. 
Spinning waveforms depend on several parameters (masses, spin components of the binary constituents, angles defining 
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the orientation of the source with respect to the observer) . For this reason it is not practical to use these waveforms 
to build template banks, but we envision their usefulness as injection waveforms in testing existing LIGO- Virgo data 
analysis pipelines, which are currently based on non-spinning templates. 

The paper is organized as follows. In sec.lTTlthe method for the analytical waveform construction and the numerical 



simulations are illustrated. In sec. Ill the results are presented, in the form of comparison between analytically and 
numerically generated waveforms. We restricted the analysis to the main quadrupolar mode {m — 2), but this setup 
can be extended to include other modes. In sec. |IV| the conclusions that can be drawn from the present work are 
reported. 

II. THE METHOD 

The goal of the present work is to produce the complete analytical waveforms generated by the coalescence of 
spinning binary systems. The numerical waveforms used to construct and calibrate our analytical model all describe 
equal mass binary systems {mi = 7712)1 with spin magnitudes |Si| = IS2I = Q.Qm\ and S2 orthogonal to the initial 
orbital angular momentum (where 81^2 denote the binary constituent spin vectors and we posit Gn = c— 1). 

The description of the dynamics adopted here models the inspiral phase via the standard TaylorT4 PN formulas, 
see [22j for definition and comparison of different PN approximants in the spin-less case. In the non-spinning case 
the waveform is fully parametrized by amplitude and the orbital phase, but in the spinning case the effect of the 
precession of the orbital plane should also be accounted for, see |23j spin and angular momentum evolution equations. 
It is convenient to define an orbital phase ij) — J ujorb dt whose evolution is given by 

dcj) v^ dv F{v) 

'di~M' ~dt~^ dE/dv ' ^' 

where M = mi + m.2 is the total mass of the binary system, F{v) and E{v) are respectively the flux emitted and 
the energy of a circular orbit with angular frequency uJorb = v^ /M . The main gravitational wave frequency few is 
related to uiorb via few = '^orb/''^- 

In order to determine the actual waveform also the spins and the orbital angular momentum have to be dynam- 
ically evolved, see e.g. |24j . In particular by parametrizing the orbital angular momentum unit vector L as 
L = (sin i cos a, sin i sin a, cos /,) it is convenient to introduce the carrier phase '^ given by 

d^ da ,„, 

— =a;„.,-cos.-. (2) 

Numerically generated waveforms are usually decomposed in spherical harmonics, in particular the five quadrupolar 
modes {I = 2) are the only non-vanishing at the lowest order in v, and the I = 2,m, = ±2 mode are the dominant 
ones. The actual shape of the / = 2, m = 2 mode, the only one which will be used here for comparison with numerical 
relativity results, is given by the following [53] 



/i^7^^ = -2'^y^[cos4(i/2)cos(2(^I' + a)) + sin4(i/2)cos(2(vI/-a)) + 0(t;)] , (3) 

where v = mim2/M^ is the symmetric mass ratio, d the source-observer distance and this formula refers to the 
inspiral waveform only. 

The functions F{v) and E{v), necessary to determine the orbital phase, are known up to 3.5PN order as far as orbital 
effects are concerned, and up to 2.5PN and 2PN level for respectively Si,2L and S1S2, SiSi, S2S2 interactions. 

According to studies in the non-spinning case [HI [25l [26] , the Taylor T4 appears to be a very good approximant up 
to a frequency ui — irfaw — 0.1/Af for the equal mass case, even though its faithfulness seems to worsen for different 
mass-ratios. 

The PN evolution is halted when LOorb reaches the value uim that is determined by comparison with numerical 
waveforms. For ujorb > <^m {^orb is monotonic here) the angular frequency is evolved according to 

^orb{t) = :; 7777;- + ^0, ^7n < UJorb and t < Uac , (4) 

1 - t/lA 

where the three unknown parameters ujo.i and Ta are fixed by requiring that LJorbjUJorb and ujorb be continuous at the 
matching point defined by Um, and trac will be defined shortly. In this phenomenological phase the amplitude of the 
waveform is evolved according to the lowest order formula eq.®, with l and a frozen to their values at the instant 



of time when ujorb — ^m (and keeping as usual v = (Alujorb)^ ) |34) . 

After this phenomenological phase, in order to smoothen the evolution of cjorb: which should settle at the ring-down 
value utrd, it turned out useful to join the time evolution of eq.Q with a smoothing function 

t ^' 



UJorbit) = UJrd - i'i^ - 7f^ j ^2 Uac < t d-Tld Uorb < ^Tnrd = 'i-^ ^rd (5) 

where tracT^2 are fixed by requiring continuity oiujorb and ujorb- For uJorb > ^rnrd = 0.8 w^d the waveforms is described 
by the damped exponential 

^2^^ = e~^''"[Acos{u:rdt) + Bsin{urdt)] , (6) 

where A and B are constant to be determined requiring the continuity of /i2.2- For each multipolar mode defined by 
a pair Z,m there is an infinity of overtones with increasing damping factor, but for our practical purposes retaining 
only the first overtone is enough. 

Note that the degree of continuity of the waveform across the phenomenological -ring-down phase transition is not 
affected by changing the value of uj^rdi but only by the number of overtones admitted to describe the ring-down phase 
in eq. (l6]), as it happens in the EOBNR construction [22]. 

The values of the ring-down frequencies and damping factors of the three lowest overtones of the / < 4 modes can 
be read from [55] as a function of the mass and spin of the final object created by the merger of the binary system. 
The final mass is determined by the algebraic sum of the constituents' masses and the negative binding energy once 
LOm is reached, and the final spin according to the phenomenological formula given in [29] . 

The damped exponential is attached when ujorb reaches 80% of the Wrd value as it turned out that this specific value 
allowed the best overlap between analytical and numerical waveform, irrespectively of the initial spins. 

The analytical waveforms just described has been quantitatively confronted with the numerically generated ones 
by computing the overlap integral 



l-hiM 



jh,{f)hl{f)df, (7) 



maximized over initial phase and time of arrival. The angular frequency LOm has been determined by a first set 
of numerically produced short waveforms (4-6 cycles long) by picking the value which allowed a maximum overlap 
integral with the phenomenological waveform, with a determination precision of ±10~''/M. 

The numerical waveforms used in the present work have been generated with MayaKranc. The grid structure for 
each run consisted of 10 levels of refinement provided by CARPET [30 , a mesh refinement package for CACTUS |31j . 
Sixth-order spatial finite differencing was used with the BSSN equations implemented with KrEoic [32] • The outer 
boundaries are located at 317M and the finest resolution is Mjll . Waveforms were extracted at 75M. A few 
waveforms were generated at resolutions of {M/64, Af/77, Af/90}, and convergence consistent with our fourth order 
code is found. The short (long) runs showed a phase error on the order of 5 • 10~^ (5 • 10~^) radians and an amplitude 
error of « 2% (w 5%). We expect similar accuracy in all runs performed. 

The numerical waveforms describe the Z = 2, to = 2 mode and consist of two sets, both with equal mass and spin 
magnitudes |Si| = IS2I = ^X>rn\. The first set consisted in 24 few-cycle-long waveforms, representing mostly the 
merger and ring-down phases of a coalescence. They have been used to fix the values of Wm for the corresponding 
values of initial spins. One of the binary constituent had initial spin 82/™^ = (—0.6,0,0) in the reference frame in 
which the initial L — (0, 0, 1). The difi'erent values of the first dimension-less spin have been obtained by rotating the 
(0, 0, 0.6) vector by 15 degrees in the x-z plane. Once determined the values of iOm for each waveform, the value of the 
matching frequency aj„i for generic spins has been determined by assuming an analytical dependence of a;„ on 81,2, 
according to 



Ma;,„ = ao -f ai(S'i^ + 5'2^) -f a2(5(S'i2 - S'2z) + a3(5i52)_L +a4(S'f_L + 52_l)-|- 
a^6{Sl^ - Sl^) -K a^{Sl, + 5|J + 07(^1.52.) + a^KSl. ~ ^L) + • ■ • 



(8) 



where the suffix _L stands for projection onto the plane perpendicular to L and higher power of the spin variables 
have been neglected. The spin components are understood in a frame where the orbital angular momentum is along 
the z-axis and their values are time-dependent, thus inducing a mild time dependence in ujm- Note that because of 
the dependence of the uJorb evolution equation on L, the spin components parallel to the orbital angular momentum 
enter already at linear level, whereas the dependence on the spin components in the plane of the orbit starts from the 



quadratic level. The Oi coefficients may depend on the symmetric mass ratio v^ but it is assumed here that they can 
be analytically expanded around their value at (5 = \/\ — Ai/ = (mi — 1712)/ {mi + 7712) = 0, according to 

a,{S) = a, + 6 4^^ + 6\l^'> + . . . . (9) 

Since the set of simulations analyzed in this work do not include the case of unequal masses, it has been possible to 
fit only the first term in the expansion ^ for all the coefficients. [3S] 

Given the specifics of the simulations we used (all having 6 = and S2 • L = initially) , it has been possible to 
determine only the coefficients oq, ai, 03, 04, ag. These values have then been tested against a second set of numerical 
waveforms, consisting of 8 long waveforms (12-15 cycles long), where all the three phases (inspiral, merger and ring- 
down) play a role. In particular in the long waveform case, even with the initial condition S2 • L = 0, the dynamical 
evolution has lead to generic S2 at the time of the onset of the phenomenological phase. Thus the determination of 
the unknown coefficients 02.7,8,... will be necessary for an accurate estimation of a;,„. Additional numerical simulations 
are required to obtain these extra information and thus improve our model. 
The result of the comparison between analytical and numerical waveform is the subject of the next section. 

III. RESULTS 

The analytical waveforms have been calibrated by comparing them with 24 short numerical simulations, with the 
result qualitatively shown in figs. l|2 and quantitatively reported in tab. IT] 



The determination of the Wm's giving the best overlap for different spin values allowed to evaluate some of the 
coefficients in the phenomenological formula (IS]), as given in tab. In] 

The coefficients 02,5,7,8 cannot be determined by the analysis performed because the terms they multiply vanish 
identically for the 24 simulations considered. The oq coefficient has been determined by comparison with a non- 
spinning waveform |36j. 

Having fixed the value of ujm with some generality, it is possible to generate analytical waveforms with any specific 
initial condition without tuning any single parameter: the value of tOm will be determined analytically via eq. (l8| with 
the unknown coefficients arbitrarily set to zero. It is then possible to generate waveforms with no tunable parameters 
for comparison with the second set of long, numerically generated waveforms, even though a wider range of initial 
conditions would be needed for a more accurate determination of u;,„ in the most generic case. Such comparison is 
summarized qualitatively in fig. [3] and quantitatively in tab. Ill 



IV. CONCLUSIONS 

We presented an analytical method to produce complete gravitational waveforms from spinning coalescing binaries. 
The free parameters of the model are the values of the orbital frequency at the transition from the inspiral to the 
phenomenological phase and from the phenomenological to the ring-down phase. They have undergone a calibration 
process by confronting the analytical waveforms with numerically generated ones and in particular the accurate 
determination of the first of such parameters {uj,n) turned to be crucial for a satisfactory waveform construction. 
For the calibration process a first set of (short) waveforms has been employed, obtaining overlap factors ranging 
between 0.97 and 0.99. Once the calibration has been obtained, the analytical waveforms, with no parameter to 
tune, had been confronted with a second set of (long) numerical waveforms and overlap factors ranging from 0.95 to 
0.99 have been obtained. Better overlap is expected once a larger number of available numerical waveforms makes it 
possible to have a more solid estimation of the matching frequency Wm for generic spins. 

The method illustrated here can be extended to produce other modes than the I = 2,m = 2 for which it has been 
tested. The m = ±2 quadrupolar modes are the only non-vanishing modes for head-on observation of a coalescing 
system, whereas to build physical waveforms for generic angle between the source and the observer all other modes 
are required, at least with I < 6. In particular, as the quality of the numerical waveforms describing such modes 
improves, it would be interesting to verify if the calibration performed here on the dominant quadrupolar mode allow 
a good description of other modes. 

The programs generating these phenomenological waveforms are written in the C language and are available from the 
LIGO Analysis Library (LAL) [33]. 
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Figure 1: Comparison between the / = 2,m = 2 mode of analytical (black, solid) and numerical (red, dashed) equal mass 
simulations: time-domain waveforms are shown in the upper frame, the lower frame shows uJorb (given respectively by its PN 
value, by eq.Q and by cord in the three different phases) vs. the instantaneous frequency of the numerical waveform (computed 
as Im(/i+ + ihx)/\h+ + ihx\, being h+ {hx ) the plus (cross) polarization of the numerical waveform). Initial spin configuration; 
Si/ttii = (0,0,0.6), S2/wi2 ~ (—0.6,0,0). Vertical dashed lines mark the onset of the phenomenological phase parametrized 
by eq. Q and the ring down phase parametrized by eq. Q. The resulting overlap integral is 0.986. 
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Table I: Values of the overlap integral between the analytical and the 24 short numerical waveforms. For reference, the values 
of iUm maximizing the overlap and of cUrd are reported. 
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Figure 2: Summary of the comparison between the I — 2,m — 2 mode of analytical (black, solid) and numerical (red, dashed) 
waveforms. Initial spin configuration: 82/77x2 = (—0.6,0,0), Si lies in the x-z plane. A rotation of Si by 15 degrees in the x-z 
plane takes from one plot to the following one, the first plot is the same as in fig. [l] Vertical dashed lines, x- and y-axis are 
the same as in fig. [l] 
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Table II: Coefficients of eq.® as determined by comparison of the analytical waveforms with 24 short numerical waveforms. 
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Figure 3: Comparison between the Z = 2,Tn = 2 mode of analytical and long numerical waveforms (mi = m2). Initial spin 
configuration: S2/m| — (—0.6,0,0), Si lies in the x-z plane. A rotation of Si by 45 degrees in the x-z plane moves the initial 
condition of one plot to the following one, the first plot refers to the initial condition Si /ml — (0, 0, 0.6). Vertical dashed lines 
mark the onset of the phenomenological phase and the ring-down phase, x-axis is time in units of M. 
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Table III: Values of the overlap integral between the analytical and the 8 long numerical waveforms. For comparison with the 
"short" case the values of ujm and ujrd are reported. 
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